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The quantum magnetic oscillations of electrical (Shubnikov de Haas effect) and thermal conduc- 
tivities are studied for graphene which represents a distinctive example of planar systems with a 
linear, Dirac-like spectrum of quasiparticle excitations. We show that if a utmost care was taken to 
separate electron and phonon contributions in the thermal conductivity, the oscillations of electron 
thermal conductivity, n(B) and Lorenz number, L(B) would be observable in the low-field (less than 
a few Teslas) regime. 
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I. INTRODUCTION 



There is a variety of condensed matters systems that in the first approximation can be regarded as planar and which 
in the low energy sector can be described by the Dirac-like form of the effective Hamiltonian. The difference in the 
behavior of the particles with the usual parabolic spectrum and the linear, Dirac-like spectrum becomes particularly 
""O ■ prominent when an external quantizing magnetic field B is applied perpendicularly to the plane. The energies of 
Landau levels for free nonrelativistic electrons are E n — ehB/(mc)(n + i), while for the "relativistic" problem 

o ■ 
o 



E n = ^Jehv F B2n/c, (1) 

co ! 

with n = 0, 1, Here e is the electron charge, m is the effective mass of carriers in the parabolic band, vf is the 

■ Fermi velocity of the system with linear dispersion, and we wrote H and c explicitly; in the following sections we set 
OO \ h = c = k.B = 1, unless stated explicitly otherwise. 

The most important qualitative difference between these two spectrums is that for the realistic values of the param- 
eters m and Vp rather weak fields B ~ 10 Tesla is sufficient to drive "relativistic" systems in the extreme quantum 
regime If causing such interesting phenomena as quantum magnetoresistance. Another qualitatively distinguishable 
feature of the Dirac-like quasiparticles is an unusual form of the semiclassical quantization condition for energy levels 
in the magnetic field S(e) = 2neB / (hc)(n + 7), where S is the cross-sectional area of the orbit in k space, n is a large 
integer (n > 0), and 7 is a constant (0 < 7 < 1). For the parabolic dispersion 7 = 1/2, which is commonly used in 
describing magnetic oscillation (MO) phenomena in metals Q, so that the corresponding Berry's phase, 7— 1/2 is zero. 
1 1 However, for the Dirac quasiparticles Berry's phase is nontrivial and 7 = 0, so that the commonly used expressions 
have to be modified accordingly. This was indeed obtained in our previous paper £| , where we have studied MO of 
the density of states (DOS), thermodynamic potential, and magnetization in QED2+1 with the spectrum QJ. We have 
also discussed the underlying condensed matter models and systems that in the low-energy approximation are reduced 
to QED2+1 form. Among them is graphite, which is probably the most promising and convenient for experimental 
investigation of the ultraquantum regime, when only a few lowest Landau levels are occupied. This semimetal was 
originally studied almost 50 years ago and there was a considerable renewal of the interest in the electronic properties 
of this material during the past decade due to the discovery of novel carbon-based materials such as carbon nanotubes 
constructed from wrapped graphene sheets [f| . While graphite itself is a three-dimensional material in which planar 
sheets of carbon atoms are stacked, graphene is an individual layer or sheet made from the carbon atoms. 

The simplified QED2+1 description is obviously appropriate for graphene. In Kish graphite the anisotropy of the 
resistivity p c along the c-axis direction and the basal-plane resistivity, pb can be p c / Pb ~ 10 2 and even as large as 
Pel Pb ~ 5 x 10 4 in highly oriented pyrolytic graphite (HOPG) @, indicating that the layers are weakly coupled. 
Therefore in the first approximation, the QED2+1 description with some limitations may also be used for graphite. 
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Indeed, early theoretical investigations of graphite show that while low- lying Landau levels correspond to Eq. QJ, 
there are other levels related to the warping of the graphite Fermi surface Q ■ Very recent measurements of de Haas 
van Alphen (dHvA) and Shubnikov de Haas (SdH) oscillations in HOPG [3] (for an earlier literature see Refs. pillflj) 
confirm that among other carriers in graphite there is a majority holes with tw-dimensional (2D) Dirac-like spectrum. 
The dependences of the thermal and Hall conductivities on the applied magnetic field in HOPG were studied in 
Ref. hj and more comprehensive data on thermal conductivity and the deviations from the Wiedemann-Franz (WF) 
law were reported in [L|. We mention aslo recent STS observations [L| of Landau levels at graphite surfaces. Finally 
we refer to monocrystalline graphitic films made by repeated peeling of small mesas of HOPG 141. These films contain 
only a few layers of graphene. The SdH oscillations are clearly observed in this material [lj| and they depend only 
on the perpendicular component of the applied magnetic field. This proves the 2D nature of the material. Moreover, 
the carrier density (and even the character of carriers, either electrons or holes) in this system is controllable by 
electric-field doping, so that there are SdH oscillations with varying applied voltage [iM- 

The purpose of the present paper is to extend the analysis of our previous paper [4J and study the MO of the 
electrical, u(B) (SdH effect) and thermal k(B) conductivities. We demonstrate that since the field dependences of 
cr{B) and k(B) at low but still finite temperatures are different, there is a violation of the WF law with observable 
oscillations of the Lorenz number 



L{B,T) = ^±-. (2) 

We begin by presenting in Sec. [H] the model Lagrangian describing Dirac quasiparticles in graphene. The general 
expressions for electrical and thermal conductivities are written in terms of the same spectral function A(B) and the 
thermal conductivity is considered also including the thermal power term. In Sec. lIIII we analytically extract magnetic 
oscillations contained in the spectral function A(B). In the Discussion, Sec. II VI our final results are summarized and 
their applicability for the graphite is considered. 



II. MODEL AND GENERAL EXPRESSIONS FOR ELECTRICAL AND THERMAL CONDUCTIVITIES 

A. Model Lagrangian 



The Lagrangian density of noninteracting quasiparticles in a single graphene sheet in the continuum limit reads 
(see e.g. 



agrangia 

IB El) 



C = r ) 



CT=±1 



V F 



H 1 (d x - ieA x ) - ij 2 (d y - ieA y ) 



*a(*,r), (3) 



where the four-component Dirac spinor ^ „ = (ipio-i tyio) is combined from two spinors tpia-, fp2tr that describe the Bloch 
states residing on the two different sublattices of the biparticle hexagonal lattice of the graphene sheet, and <j = ±1 
is the spin. In Eq. © ^ a = is the Dirac conjugated spinor and 4 x 4j matrices are either (0-3, io"3, —iai) <8> 03 

[ls| or their unitary- equivalent representation can be taken from Ref. 0, 0] . 

Note that there is no principal difference between between two-band models for electron and holes discussed in 
Refs. [53, El and a model with Dirac fermions, where these electron and holes with identical velocities are built in 
the formalism. There are, however, some cases like a double-resonant Raman scattering in graphite 22], where the 
asymmetry between the bonding, E = Vp\k\ and antibonding E = — v' F \k\ bands in graphite is essential, so that 
multiband models are more suitable. 

Since the terms with d x . v in Eq. originate from the usual kinetic term of the tight-binding Hamiltonian, vector 
potential A is inserted in the Lagrangian using a minimal coupling prescription. The vector potential for the 
external magnetic field B perpendicular to the plane is taken in the symmetric gauge 

B B \ 

7X2,— x x ). (4) 



2 2 

Using for the value of the nearest-neighbor hopping matrix element of graphite t ~ 2.3eV, we obtain that the Fermi 
velocity is vf ~ 7.4 x 10 5 m/s, and accordingly one can estimate from Eq. that E\ ~ 300-^ • \J B [Tesla] . 

Since the Lagrangian © originates from nonrelativistic many-body theory, the interaction of the spin degree of 
freedom with magnetic field 

L B =MB£^<7f <r (t,r) 7 °* CT (t ) r) (5) 
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has to be explicitly included by considering spin splitting fx n = /j — ct/j-bB 2] of the chemical potential /i, where 
fiB = eh/(2mc) is the Bohr magneton. Note that the number of spin components can be regarded an additional 
adjustable flavor index of fermions a = 1, . . . , N and N = 2 corresponds to the physical case. The magnitude of the 
Zeeman term depends on the ratio fj-s/ks — 0.67K • Tesla -1 . This term, in fact, has the same magnitude as the 
distance between Landau levels in the nonrelativistic problem. Although we will include this term for completeness 
in the analytical expressions, in the numerical calculations it can be safely neglected because it is much smaller than 
estimated above E\. 

To make the treatment more general, we also include a mass (gap) term 

C A = ]T <M*,r)A<Mt,r), (6) 

<J=±\ 

in the Lagrangian J2J). For example, it is well known that an external magnetic field is a strong catalyst in generating 
such a gap for Dirac fermions (the phenomenon of magnetic catalysis) 23]. Usually the opening of the gap marks 
an important transition which occurs in the system. In particular, in the case of pyrolytic graphite a poor screening 
of the Coulomb interaction may lead to excitonic instability, resulting in the opening of the gap in the electronic 
spectrum and manifesting itself through the onset of an insulating charge density wave (see e.g. [HHIi^ . l24 | ) . 



B. Electrical conductivity 



The dc conductivity tensor can be found using Kubo formula [2£ 



hm J — , (7) 



where 11^ (fi) is the retarded current-current correlation function (see e.g. [l9ll26j 'l 

p . e 2 vl v-^ f°° , , tanh[(w 2 -<u CT )/2T] - tanh[(wi - U a )/2T\ f d 2 k rA , 

J 2 ^ J_ QO uji - lu 2 + £1 + tO J (2tt) z 

Here A(uj, k) is the spectral function associated with the translationary-invariant part of the Green's function of Dirac 
quasiparticles in an external magnetic field given by Eqs. (3.3) and (3.6) of Ref. || (see also Refs. EHHHl). Then 
for the diagonal conductivity a = a xx — a yy we have 

a = ne 2 v 2 F J2 J J <M-nJp(w - [A{uj, k) 7l A(u, k) 7l ] , (9) 

a — oo 

where — n' F (uj — //) = (1/4T) cosh -2 [(a; — fj,)/2T] is the derivative of the Fermi distribution. Further details of 
calculation of tr and the momentum integral in Eq. © were considered in Refs. [lH, l26l l27j . so here we write down a 
rather simple final expression for the electrical conductivity in an external magnetic field 

oo 

a = e 2 Y [ % A(w,B,T,A), (10) 

^ I 4Tcosh 2 ^^ V ' ' ; ' 



where the function 



At R r a n - 1 P2 f 2 , (^ 2 + A 2 + T 2 )(eB) 2 - 2u\u 2 - A 2 + T 2 )eB 

a, , aj - ^ (gB)2 + (2 ^ r)2 ^ + (w3 _ A2 _ r2)2 + 4w2p2 

w (^ 2 -A 2 + r 2 ) /A 2 + r 2 -^ 2 -2^r 

f (, 2^ 



(11) 



Here in order to consider the MO for a more realistic case, we introduced the effect of quasiparticle scattering by 
making 5-like quasiparticle peaks associated with the Landau levels Lorentzians with a constant energy-independent 
width r (see details in Ref. 4]). This approximation still allows us to derive a rather simple analytical expression 
(|lip. where ip is the digamma function, which eventually results in a Dingle factor in the expression for the amplitude 
of MO. 
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Note that in general one should consider dressed fermion propagators that include the self-energy S(w) due to the 
scattering from impurities. Up to now the problem of scattering from impurities in the presence of a magnetic field 
does not have yet a satisfactory solution. Therefore, here we have chosen the case of constant width r = T(lu = 0) = 
-Im£ R (u> = 0) = l/(2r), r being the mean free time of quasiparticles. 

Such a Lorentzian broadening of Landau levels with a constant T was found to be a rather good approximation 
valid in not very strong magnetic fields 0, • Definitely, the treatment of disorder in the presence of the magnetic 
field in such a simplified manner should be considered as only a first step until further progress in this problem is 
achieved (in connection with this, see, Refs. |2g.l30|'). 



C. Thermal conductivity 

The longitudinal thermal conductivity can also be calculated using a thermal Kubo formula [2£ 

k(B, T) 



T 



_L_ Imng g (n + iO) 

t 2 n 



i lim [imn| c (o + i o)] 2 



T 2 a n^o 



ft 2 



(12) 



Here H^ E (fl) is the retarded longitudinal energy current-current correlation function and Il|| c (f2) is the retarded 
longitudinal correlation function of energy current with electrical current and the expression for the energy current 
operator is given in Refs. [2(| |53 ■ 

We note in passing that the second term of Eq. (|12fl is related to the thermal power 



S 



1 , lmn§ r (Q + iO) 

— lim z±i ±. 

Tn^o Im±T«(ft-M0) 



(13) 



where II"(f2 + iO) = H xx (n + iO) = Tly y (fi + iO). The presence of the thermal power term in Eq. I|12fl ensures that 
the energy current is evaluated under the condition of vanishing electrical current |25| . 
Similarly to the above-derived electrical conductivity we finally arrive at 



<B,T) 
T 



E 



dijj 



a(B,T) 



U - Ho 

T 



E 



i 



4Tcosh 2 ^ 



——A(u),B,T,A) 



T 4Tcosh 2 ^ 



—A(u,B,T,A) 



(14) 



where A is the same function l|ll|) as for the electrical conductivity. This function contains all information about the 
field dependence of the transport properties of the systems with a linear dispersion law, including the MO. While the 
representation (jll|) can already be used for numerical calculations, for analytical work it is useful to extract explicitly 
the MO that are contained in the digamma function ip when the real part of its argument becomes negative. 



III. ANALYTICAL CONSIDERATION OF OSCILLATIONS 



A. Extracting oscillations from A using ?/)-function properties 



The oscillations of A in 1/ B can be extracted using the relationship for ip function 

ip(—z) — i/j(z) H h 7T cot irz, 



which [see also Eq. (4.18) of Ref. jj]] results in the expression 
'A 2 - (e + iT) 2 



4> 



= Re?/> 



2eB J r V 2eB 
e 2_ A 2 -r 2 | -2ier 



/ - [ A ' ^P^- ) W - a 2 - r 2 ) + 9(A 2 + r 2 - e 2 )] 



2eB 



- i sgn(e 2 - A 2 - T 2 )lmip N 



-A 2 -r 2 |-2zer 



2eB 



6(e 2 - A 2 - T 2 



2eB 



A 2 - T 2 + 2ieT 



7T COt 7T- 



e 2 - £^T^+2ieT 
2^B 



(15) 



(16) 
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Taking the imaginary part of the last equation, we obtain 
.2 / , . i „-n\2 



ImV> ( 7—^- — - J = -sgn(w 2 - A 2 - r 2 )Im^ 

- e(u? - a 2 - r 2 



AeBuoT 



- A 2 - T 2 1 - 2iujT 
2e~B 

sinh(27rcjr/ eJ B) 



(17) 



> 2 - A 2 - T 2 ) 2 + 4j 2 r 2 cosh(27rwr/e-B) — cos[it(uj 2 — A 2 — T 2 )/eB] 
Then substituting Eq. (|T7Jl in Eq. ljTT|) we obtain 

1 T 2 f_ , (eB) 2 /2 (eB) 2 /2 2w 2 (w 2 - A 2 + r 2 )e£? 



^,B,r,A)-^ 2 (eB)2 + (2 ^ r)2 

w(u; 2 - A 2 +T 2 ) r 



2u/ 



(^ + A) 2 + r 2 (w-A) 2 + r 2 (w 2 - A 2 -r 2 ) 2 + 4w 2 r 2 
- a 2 - r 2 1 - 2i^r\ 



2eB 



-6>(w 2 



r 2 



sgn(u> 2 - A 2 - r 2 )Im^ 
AeBuY 

{lo 2 -A 2 - T 2 ) 2 + 4w 2 r 2 + 7r cosh(2 7 ra;r/eS) - cos[7r(w 2 - A 2 - r 2 )/eS] 



(18) 



sinh(27rwr/eS) 



where the oscillations are contained in the last term of Eq. (|18f) . Note that the real part of the argument of if> function 
in l|18|) is already positive and the signs before ip in l|ll|) and (|18|l are different. 
For u! 2 > A 2 + r 2 using the relationship 



- (a — ib) 



Re 



cos b — e a 



l-e-( Q - ib ) 2(cosha - cos 6)' 



one can expand 



Ao SC (uj,B,T,A) =- 



{eB) 2 + (2ujT) 2 

irk(Lu 2 - A 2 - T 2 ) / 2,T/,-U'|r 
x > cos exp 



k=l 



eB 



eB 



(19) 



-, sinh(2^nr/e5) = 1 + 2 V cos[7rfc(^ 2 - A 2 - T 2 )/eB] exp(-2^fc| W |r/eS) (20) 

cosh(27ro;r/eS) - cos[7r(w 2 - A 2 - T 2 )/eS] ^ V ; 

and finally arrive at the expression for oscillatory part of A 

2 ujY{uj 2 - A 2 + T 2 )6(u; 2 -A 2 - T 2 ) 



(21) 



1. Low-field non oscillatory limit 
Eq. (|18fl can be simplified in the low field limit: 

1 - 

+6{lu 2 - A 2 - r 2 



A{^B,Y,A) = 1 1 



u 2 -A 2 +T 2 , o A , 2 , 2Mr 

sign(w 2 — A — T ) arctan — 



2|w|r 

tt(lu 2 -a 2 + r 2 ) 



-A 2 -r 2 | 

sinh(27rwr/eB) 



2ujT cosh(27rwr/ eJ B) - cos[tt(w 2 - A 2 - T 2 )/eB] ' 

where we kept B only in the oscillatory part of A. For uo 2 < A 2 + T 2 after using the relationship 

x > 0, 



7T 1 

arctan x = arctan — , 

2 x 



the last equation reduces to [see also Eq. (4.16) of Ref. pfij ] 



A(oj,B = 0,r,A) 



2tt 2 



-A 2 + r 2 /tt A 2 + r 2 -^ 2 

— ; — arctan ; — r- 

2Ur v 2 2 \u r 



(22) 



(23) 



(24) 
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B. Oscillating parts of electrical and thermal conductivities 

Substituting Eq. (|21|l in (|10fl one can obtain the expression for oscillating part of conductivity 

_4e 2 Hr (m 2 - a 2 + r 2 )#(/i 2 - a 2 - r 2 ) 



x cos 

fc=l 



{eB) 2 + (2^F) 2 
7rfc(/i 2 — A 2 — r 2 ) 



eB 



(25) 



i?T(fc, /i)i?D(fc, m), 



where we introduced the temperature amplitude factor 



R T (k,fj,)=R T (t k ) = - 



sinh tfe ' 



2ir 2 kTn 
eB 



(R T (0) = 1), 



and the Dingle factor 



-Ro(fc, /i) = exp 



27r%|r 



eB 



(26) 



(27) 



Deriving Eq. i|25[l we made the following simplifying assumptions: (i) Spin splitting is not included, and (ii) the low 
temperature T — > limit is considered. Thus after making a shift lj — > + and changing the variable w — > 2Tw, 
and keeping only the linear in T terms in the oscillating part of the integrand, we used the integral 



dx 



cos free irb/2 
cosh 2 x sinh7r&/2 



(28) 



to obtain the temperature amplitude factor l|26|) . It is essential that in contrast to the the Dingle and temperature 
factors for nonrelativisitic spectrum, both (|26[) and l|27|l factors for the relativistic spectrum contain chemical potential 
fi (see also Ref. Q). The distinctive concentration dependence of Rt should be observed experimentally. 
Similarly to Eq. 1251) one arrives at the expression for the oscillating part of thermal conductivity l|14l) : 



4ir|At|r {n 2 - A 2 + r 2 )^ 2 - A 2 - T 2 ) 



T 
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x > cos 

k=l 



(eB) 2 + (2//F) 2 
nk(^ 2 - A 2 - T 2 



eB 



(29) 



R^(k,n)RD(k,fj,), 



where 



R'x(tk) 



12 d 2 R T (t k ) 



G 



sinhtjr. 



cotht fe - — 

2 sinh t k 



[Rt(0) = 1] 



(30) 



is the temperature amplitude factor for thermal conductivity obtained using the second derivative of the integral i|28[) . 
The temperature amplitude factors Rr(t) and Rt^(t) are shown in Fig. 2] Interestingly the dependence R^{t) is 
non-monotonic and R^(t) even changes its sign. 



IV. DISCUSSION 



Based on Eqs. (fTUf> and i|14fl. in Figs. [21 [SI we compute the field dependences of a(B,T)/a(B = 0, T), and 
k(B,T)/ k(B = 0,T) and the normalized Lorenz number L(B,T)/Lq, where Lq = tt 2 k 2 B / (3e 2 ) is the Sommerfeld's 
value for the Lorenz ratio. To simplify our consideration we set A = 0, but it may be necessary and quite interesting 
to consider the influence of A on MO; see Ref. 0. As mentioned above Eq. J()J|, we also do not include spin splitting. 
Analyzing these figures one can conclude the following. 

1. As usually, the conditions favorable for the magnetic oscillations are T,T <C w^, where ljl ~ E\ is the dis- 
tance between Landau levels. This regime is different from the regime T <C ujl < T of an unconventional 
magnetotransport [2lJ when MO's are still not resolved due to the thermal smearing of Landau levels. 
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2. As the field increases from B = both electrical and thermal conductivities decrease rapidly and start to 
oscillate when Landau levels cross the Fermi surface. The specific of graphene is that for the realistic values 
of the parameters there are only a few Landau levels below the Fermi surface. As the field B increases, these 
levels quickly cross the Fermi surface, and as one can see from Fig. |3 for B > 4Tesla the lowest Landau level 
Ei is already above the Fermi surface so that the MO's disappear and the system can, in principle, enter in the 
quantum Hall effect regime. Thus MO's in graphene with eB < /it 2 are quite different from the conventional 
MO's when there are so many Landau levels below the Fermi surface that u>l /i and when a new level crosses 
the Fermi surface the system returns practically in the same state as before. 

For Dirac fermions in order to have at least one oscillation the inequality [see Eq. (8.20) of jj|] /it 2 — A 2 > 2eB 
have to be satisfied, and analyzing the experimental data of Refs. 0, ITlT ] there is a temptation to distinguish 
there the regimes of conventional MO's and quantum Hall effect. However, the real situation is more involved 
due to the difference between the model for graphene with linearized spectrum and measured properties of bulk 
graphite (see also Ref. 01). Nevertheless, the observation of plateaulike features in the Hall resistivity for 
B > 2Tesla |ll| suggests that for these fields the system is already in the quantum Hall effect regime. 



3. Comparing Figs. and [31 one can see that cr(B) and n{B) do not oscillate in phase. In particular, we observe 
that at T = 28K the oscillations of a(B) are practically invisible, while one can still notice some oscillations of 



To look closer at the difference between a and n in Fig. 01 we also plotted the thermal conductivity kq(B,T) 
calculated without the second term of Eq. (|12|l . In normal metals this second term is considered to be unimpor- 
tant because usually it is ~ T 2 //u, 2 times less than the first term of Eq. l|12f) and because the WF law is always 
considered in the limit T <C fJ. this term is usually neglected. 

4. By comparing Fig. |3 and Fig. 0] it is easy to see that each peak of the electrical conductivity is accompanied 
by two satellite peaks of the thermal conductivity. The dip between these two peaks in k(B) coincides with 
the peak of <j(B). The origin of these satellite peaks is related to the fact that the expression for the thermal 
conductivity l|12f> contains the factor g(u>) = —n' F (uj — fJ,)(u> — /i) 2 /T 2 and thus measures A(u>, B,T) below and 
above the Fermi energy, while the electrical conductivity probes A[ui, B, T) at the Fermi energy, because it 
contains just the factor n' F (u>) (see j3]| and Figs. 12 and 13 of Ref. |2q|L 

5. Let us now compare Figs. |3 and 0] At T = 3K each double-peak structure observed in Fig. 01 is replaced by a 
single broader peak in Fig.|3 This reflects the fact that the full expression (|12f> for k(B) ~ kq(B) — a(B)/a(B), 
so that if the coefficient a(B) is large enough, the peaks seen in <r(B) also produce an increase of k(B), so the 
dip between peaks in k(B) is filled in and we observe a single broader peak. 

For higher temperatures the role played by second term of Eq. (|12fl further increases and k(B) and Kq(B) behave 
quite differently. Finally we observe the above-mentioned picture when at T = 28K the oscillations of cr(B) are 
damped, but one can still see some oscillations of k(B), but the positions of the peaks do not coincide with the 
positions of the peaks in <j(B) and k(B) observed at lower temperatures. 

6. The behavior of Lorenz number is shown in Fig. [S] Since in the chosen temperature interval the inequality 
T -C /i is well justified, we observe that the WF law is maintained in zero field (see Ref. 26]). In nonzero field 
we observe violations of the WF law that become more strong as the temperature increases. At T = 3K the 
behavior of L(B) is similar to the behavior of kq(B) and two satellite peaks in L(B) are related to the broad 
peak in k(B), while the dip between these peaks is caused by the peak of a{B). As the temperature increases 
a more complicated behavior of k(B) results in the large-amplitude oscillations of L(B). The positions of the 
peaks are not related to the positions of low-T SdH oscillations. In the high-field regime B > lOTcsla there is a 
tendency to the restoration of the WF law. 

The above-presented picture is already quite complicated due to the interplay between the first and second terms 
of Eq. i|12|) , so that in general there are no correlations between the low temperature SdH oscillations and oscillations 
in k(B) and L(B) seen at higher temperatures. Further complications can be caused by the fact that the impurity 
scattering rate T which we assumed to be field and temperature independent may, in fact, depend on both B and T 



the results for electrical conductivity shown in Fig. |21 one should bear in mind the difference between graphene and 
graphite mentioned in item 2. In Refjg the measurements were done at T = 2K and the oscillations are clearly seen 
and some of them, as stated in Ref. y|, are related to the q uasip articles with a linear dispersion. It is likely that for 
higher temperatures the oscillations of a(B) are suppressed [ll| and we also see that at T = 28K the oscillations of 
<t(B) are practically smeared out. 



k(B). 
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A comparison with experimental results for thermal conductivity is more complicated because the measured thermal 
conductivity K exp (B, T) = n(B, T)+k p ^(B, T) besides the electron contribution, k(B, T) contains also the contribution 
from phonons, K p b(T), which is assumed to be field independent. Accordingly, we can relate theoretically calculated 
k(B,T)/ k(B = 0,T) shown in Fig.[31with the experimentally accessible ratio K cxp (B,T)/K exp (B = 0,T) via 

K cxp (B,T) K (B,T)/k(B = 0,T)-1 K (B,T)/k(B = 0,T)-1 

K exp (B = 0,T) + Kph (T)/n(B = 0,T) + l ~ + (T/T ) 2 + 1 ' [ ' 

The last equality is written using that the ratio Knh(T)/ 'k(B = 0,T) can be determined from the fact that at 
T = T = 1.5K the values n ph (T ) ~ n(B = 0,T ) so that assuming n ph (T) ~ T 3 and n(B = 0,T) - T we 
estimate K ph (T)/n{B = 0,T) ~ (T/T ) 2 . 

The ratio K exp (B,T)/ K exp (B = 0, T) is plotted in Fig. El As one can see the MO's of n(B,T) are masked by 
the phonon contribution in n exp (B ,T) and only at rather low temperatures is there a possibility to observe them. 
Moreover, since the amplitude of MO is higher for lower fields (see also Figs. 1214^ . the MO's are more easily observable 
in the low-field regime. 

Taking these facts into account let us discuss the oscillations of thermal conductivity that were experimentally 
observed in Refs. 0, 0] for B > 2Tesla. Since these oscillations are clearly seen only in the high field regime and 
their amplitude increases as the field grows, their origin is unlikely related to the conventional SdH-like MO's and 
seems to be more associated to the quantum Hall effect. It would be interesting to check whether the conventional 
MO's of the electronic thermal conductivity n(B, T) predicted in this paper can be also observed in graphite when 
the measurements will be done at sufficiently low temperatures and/or the phonon contribution is subtracted. 
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FIG. 1: (Color online) The temperature amplitude factors Rr(t) and Rr{t) given by Eqs. JgjJ and lEiUt . 
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FIG. 2: (Color online) The normalized conductivity a(B ,T)/o(B = 0, T) as a function of field B for three different values of 
temperature T for fi = 600K and T = IK. We use eB -> (4.5 x 10 4 K 2 )B(Tesla). 
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FIG. 3: (Color online) The normalized thermal conductivity k(B,T)/k,{B — 0, T) as a function of field B for three different 
values of temperature T for /j, = 600K and T = IK. We use eB -> (4.5 x 10 4 K 2 )B(Tesla). 




FIG. 4: (Color online) The normalized thermal conductivity kq(B,T)/k(B = 0,T) [calculated without the second term of 
Eq. (II 2| that originates from the condition of absence electrical current in the system) as a function of field, B, for three 
different values of temperature T for fi = 600K and T = IK. We use eB -> (4.5 x 10 4 K 2 )B(Tesla). 
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FIG. 5: (Color online) The normalized Lorenz number L{B,T)ILq as a function of field B for three different values of 
temperature T for (i = 600K and V = IK. We use eB -> (4.5 x 10 4 K 2 )B(Tesla). 




FIG. 6: (Color online) The normalized thermal conductivity K exp (B,T)/K exp (B = 0, T) [see Eq. (I3H 1 as a function of field B 
for three different values of temperature T for (i = 600K and T = IK. We use eB (4.5 x 10 4 K 2 )B(Tesla) and T = 1.5K. 



